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Abstract Quantum discord Q is a function of density matrix elements. The 
domain of such a function in the case of two-qubit system with X density 
matrix may consist of three subdomains at most: two ones where the quantum 
discord is expressed in closed analytical forms (Q ^/2 and Qq) and an interme¬ 
diate subdomain for which, to extract the quantum discord Qg, it is required 
to solve in general numerically a one-dimensional minimization problem to 
find the optimal measurement angle 9 € (0,7r/2). Hence the quantum discord 
is given by a piecewise-analytic-numerical formula Q = min{Q^/ 2 : Qs: Qo}- 
Equations for determining the boundaries between these subdomains are ob¬ 
tained. The boundaries consist of bifurcation points. The Qg subdomains are 
discovered in the generalized Horodecki states, in the dynamical phase flip 
channel model, in the anisotropic spin systems at thermal equilibrium, in the 
heteronuclear dimers in an external magnetic field. We found that transitions 
between Qg subdomain and (57r/2 and Qq ones occur suddenly but contin¬ 
uously and smoothly, i.e., nonanalyticity is hidden and can be observed in 
higher derivatives of discord function. 
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1 Introduction 

At present we have a situation where further miniaturization of electronics 
will inevitably lead to molecular size components. Designing such components 
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requires application of the laws of quantum mechanics. This is expected to lead 
to a technological breakthrough which will be achieved through employing the 
holy of holies of the quantum theory — so-called quantum correlations. 

Initially, the entanglement has been considered as a quantum correlation 
m- Quantum entanglement is able to bind different parts of systems, even 
in the case when there is no interaction between those parts (the Einstein- 
Podolsky-Rosen effect). Quantum entanglement exists only in nonseparable 
states of bi- and multipartite systems. However, it appears in the last years 
that there are quantum correlations more general and more fundamental than 
entanglement. In particular, they can be present in certain separable states, 
i.e., when the quantum entanglement is absent. As a measure of total purely 
quantum correlations in bipartite systems, the quantum discord is employed 
now [31I311S] . The basis for the discord conception is the idea of measurements 
performed on a system and maximum amount of classical information being 
extracted with their help. 

Due to the fact that it is necessary to solve the optimization problem, 
the evaluation of quantum correlations, especially discord, is extremely hard 
[^. If for the two-qubit systems the quantum entanglement of formation has 
been obtained for the arbitrary density matrices [7], the analytical formulas 
for the quantum discord were proposed for X states [SlI^ fTOlITTlIT^fT^ . In an 
X matrix, nonzero entries may belong only to the main diagonal and anti¬ 
diagonal [mm- Notice that the sum and product of X matrices is again the 
X matrix (i.e., a set of X matrices is algebraically closed). 

However, it was found later that the formulas mMum are incorrect in 
general. The reason is that the authors mMsm believed (and it was their 
error) that the optimal measurements are achieved only in the limiting points, 
i.e., at the angles 0 = 0 or 7r/2 (see below). But on the explicit examples 
1171118] of X density matrices it was proved that the optimal measurements can 
take place at the intermediate angles in the interval (0,7r/2). Unfortunately, 
these examples with density matrices are specific and do not clarify the general 
situation. 

In the present paper we show that the domain of intermediate optimal 
angles can arise in the vicinity of transition from the domain with optimal 
measurement angle 0 = 7r/2 to the domain with optimal angle 0 = 0 (or in¬ 
versely). We derive the equations for the boundaries between these domains 
and investigate their solutions for different models. In particular, the bound¬ 
aries can coincide or be absent at all and then the quantum discord is given 
in the total domain of definition by closed analytical formulas. 

In the following sections, the general seven-parameters X density matrix is 
reduced to the five-parameter form by using local unitary transformations, the 
existence of intermediate subdomains with the optimal anglers 0 ^ 0,7r/2 is 
proved, the equations for boundaries between different subdomains are derived 
and then applied to various physical systems. Finally, a brief conclusion is 
given. 
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2 Real non-negative form for the X density matrices. Domain of 
definition for the matrix elements 


In the most general case, the X density matrix of two-qubit {A and R or 1 
and 2) system is given as 


PAB = 


fa 0 0 Ml -I- iu2 \ 

0 b vi + iv2 0 

0 vi — iv 2 c 0 

y Ml — iu2 0 0 ^ / 


( 1 ) 


where a-|-&-|-c-|-d = 1. This matrix contains seven real parameters which, 
due to the non-negativity definition of any density operator, must satisfy the 
inequalities 

a, 6, c, d>0, ad>Ui+U 2 , bc>vl+V 2 - (2) 

The quantum entanglement and quantum discord are invariant under the 
local unitary transformations of density matrices Owing to this 

property, one can with the help of transformation 


U = 


^-iif2cr^l2 


( 3 ) 


reduce the seven-parameters density matrix o to the real five-parameters X 
form [TSlfroil^l^ . UpabU'^ -a pab- This provides with the angles 


’ 2 Ml Ml 

After this, the density matrix m takes the form 


( 4 ) 


PAB = 


^ a 0 0 M^ 

0 6 M 0 

0 M c 0 

^M 0 0 dJ 


where 


M = Ml cos(arctan —) + U 2 sin(arctan —), 

Ml Ml 

M = Ml cos(arctan —) -I- M 2 sin(arctan —). 

Ml Ml 


( 5 ) 

( 6 ) 
( 7 ) 


Moreover, with the help of local rotations again around the z axis, it is not 
difficult to obtain also the non-negative off-diagonal elements of the X matrix 
Indeed, the local unitary transformation 


/ * 


t/i = e“ 




( 8 ) 


V -V 
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changes the sign of the off-diagonal matrix element u: 


UiPabUi = 


^ a 0 0 
0 b V 0 
One 0 
\-u00 d) 


Similarly, the local transformation 


U 2 = e* 


i 

—i 

\ 1 / 


selectively acts on the sign of v: 


U2PABU2 


^ a 0 0 u\ 

0 b —V 0 
0 —V c 0 

0 0 dJ 


(9) 


( 10 ) 


( 11 ) 


Thus, after transformation of X matrix to the real form we may simply enclose 
the off-diagonal elements in the modul symbols: 


PAB 


^ a 0 0 |u|\ 

0 & |n| 0 

0 |n| c 0 

\|m| 0 0 d J 


( 12 ) 


This operation does not influence on the value of quantum correlations in the 
system. 

Thus, one can now consider the density matrices with restrictions 


a,b,c,d,u,v > 0, a + b + c + d=l, ad>u^, bc>v^. (13) 


These relations define the domain V oi X density matrix in the space of its 
entries. 

We can rewrite the density matrix in the equivalent form 


PAB 


^1 + Sl -|-S2-|- 


C 3 


Cl - C2 


where 


Cl - C2 \ 

Si — S2 — C3 Cl -I- C2 
Cl -I- C2 1 — Si -I- S2 — C3 

1 - Si - S2 -I- C3 / 
(14) 


si = a-|-6—c — d, S 2 = a — b + c — d^ 

Cl = 2 ( 1 ;-I-u), C 2 = 2(u — m), C 3 = a — b — c + d. (15) 
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Decomposition of this matrix on the Pauli matrices ctq, (a = x, y, z) leads to 
its Bloch form 

PAB = (8) 1 + S2I 0 cr^ + CiCTa; ®CFx+ C 2 (Jy 0 (Ty + c^az 0 az). (16) 

The expansion coefficients are the unary and binary correlation functions and 
therefore five parameters of density matrix are expressed through the five 
different correlators, 

Si = (ctI) = Tr(p^scr2 0 1), S2 = (cr^) = Tt{pab^ 0 dz), 

Cl = (alal) = TT{pABfTx 0 <7x), C2 = = TT{pABcry 0 (Xy), ( 17 ) 

C 3 = (cr^cr^) = Tr{pABCrz 0 cTz). 


It is clear that 

- 1 < Si,S2,Cl,C2,C3 < 1. (18) 

The domain of definition, V, in the space (si, S 2 , Ci, C 2 , C 3 ) is formed, ac¬ 
cording to Eqs. m and m, by conditions (see also mm) 

(1 - Csf > (si - 52 )^ + (ci -b C 2 )^, (1 -b Csf > (si -b 52 )^ + (ci - C 2 )^. (19) 

The solid 27 is finite and lies in the five-dimensional hypercube (1181) . Numerical 
calculations show that the volume of 27 is 8 % of the hypercube one. 

The domain 27 is bounded by two quadratic hypersurfaces 

(Sl - 82 )^ + (Ci -b C2)^ - (C3 - 1)2 = 0 (20) 

and 

(si -b 32 )^ -b (ci — € 2 )^ — (c 3 -b 1)^ = 0. (21) 

Rotation by the angle 7 r /4 around the C 3 axis transforms these hyperquadrics 
to the forms 

(S 2 )" + (c'l)" - = 0 (22) 

and 

(,;)2 +(4)2-(£1^ = 0, (23) 

where 

s'i_2 = (±si + S2)/V2, c'l 2 = (±ci -b C2)/\/2. (24) 

Thus, the five-dimensional domain 27 results from an intersection of two conic 
hypercylinders ( 1 ^ and ( 1 ^ . 
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3 Three alternatives for the quantum discord 

As mentioned above, the measurement operations lie in the ground of discord 
notion. Following the founders of discord conception and their adher¬ 

ents EiiinniiiiKi] we will consider here the one-dimensional projective mea¬ 
surements (i. e., the von Neumann measurements). Such measurements can be 
reduced to projections which are characterized by the polar (0) and azimuthal 
(</>) angles relative to the z axis. It is important that in the case of real X 
density matrix with an additional condition uv > 0 the optimal measurements 
are achieved by cos 2(j) = 1 [iHiiin]. Since the sign of off-diagonal elements are 
changed by the local unitary transformations, we can always satisfy the above 
condition. 

In the general case, quantum discord depends on subsystem {A or B) 
where the measurements are performed. For definiteness and without loss of 
generality, let the measured subsystem be B. Then the quantum discord is 
given as P1H1I3] 

Q = S{pb) - S{pab) + min Scond{d), (25) 

u 

where pB = Tr^PAS is the reduced density matrix and S{p) = —Trplnp is 
the von Neumann entropy for the corresponding state p. (Here the entropy 
is in nats; to transform it, e.g., in bits, one should divide it by In 2.) Simple 
calculations with ([S]) lead to 

S{pb) = —(a + c) ln(a + c) — (b + d) ln(6 -|- d), (26) 

S{pab) = S, where 

a + d + J{a — dY + a + d + J[a — dY + Av? 

<5 =--In--- 

2 2 

a + d — \/ia — dY + a + d — \/ia — dY + 

--In-- 

2 2 

b + c + ^/(b — cY + , b + c + ^/(b — cY + 

--In-- 

2 2 

b + c — \/(b — cY + , b + c — \/(b — cY + 

--In--• 

2 2 

The quantum conditional entropy of subsystem A is given as m 

4 

Scond(0) = Ai In Ai -I- A 2 In A 2 - ^ A* In Ai, 

i=l 

Ai ,2 = i[l ± (a — 6 -I- c — d) COS0], 

Ai ,2 = ^ [1 + (a — 5 + c — d) cos 9 

± {[a-l-6 —c—d-|-(a — 6 — c-l-d) cos0]^ + sin^ (30) 


(27) 

(28) 

(29) 


where 
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Fig. 1 A fragment of phase diagram with three possible subdomains for the A-state quan¬ 
tum discord 


A 3^4 = -{1 — {a — b + c — d) cos 9 

± {[a + b — c — d — {a — b — c + d) cos0]^ + sin^ (31) 

w = |m| + |t;|. (32) 

Thus, Scond depends in fact on four parameters because u and v enter via the 
combination (I321). The conditional entropy Scond{9) is a differentiable function 
of its argument 9. 

Expressions (IMl)-® allow to define the measurement-dependent discord 
as [U 

Q{9) = S{pb) — S{pab) + Scond{9), (33) 

where 9 £ [0,7r/2]. It is obvious that the absolute minimum of this discord can 
be either on the bounds {9 = 0,7r/2) or at the intermediate point 9 € (0,7r/2). 
As a result, there is a choice from three possibilities for the quantum discord 

Q = min{Qo,Qe,Qf }• (34) 

This equation generalizes the one proposed earlier for the quantum discord [51 

iiiniiiiiiii] 

Q = min{Qo,Qf }, (35) 

i.e., it was assumed that the optimal observable can be either cr^ or Gx- In 
Fig. d] we schematically illustrate the parameter domain of a system with 
three possible subdomains for the discord. 

From Eqs. (uni)-®, we have for the discord branch Qo = Q(0): 


Qo = —S — a In a — blab — clnc — dhid. 


(36) 
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For 9 = 'K12 we obtain 

= —S — (a + c) ln(a + c) — {b + d) ln(6 + d) 

I + (a + b — c — dV + Aw'^ , l + -v/(a + 6—c — d)^ + 4w^ 

--In-- 

2 2 

1 — x/(a + 6 — c — + 4w^ , 1 — \/(a + b — c — dy + 4w^ . . 

- 2 -- 2 -• 

Thus, the branches Qg and QTr /2 are expressed analytically, and the branch 
Qe = niineg(o_,n-/ 2 ) <3(6*), if the intermediate minimum exists, should be found 
from the numerical solution of one-dimensional minimization problem or from 
the transcendental equation 

SLnA(^) = 0 . ( 38 ) 


The first derivative of conditional entropy with respect to 0 is equal to 

4 

SLnAff) = y (1 + InTi) + (1 + InT^) - ^ A'(l + In A.) (39) 

i=l 


with 


yl^2 = ~ ^ ^ ~ 9 , 


(40) 


^i,2 - 4 

j_ [a + b 


— {a — b + c — d) sin 9 

— c—d+{a — b — c + d) cos0][—(a — 5 — c -I- d) sin 0] -I- 2w‘^ sin 20 
^J[a + b—c — d+{a — b — c + d) cos 0]^ -I- sin^ 9 


(41) 


^ 3.4 - 4 

j_ [a + b 


{a — b + c — d) sin 0 

— c—d — {a — b — c + d) cos 0] (a — 5 — c -I- d) sin 0 -|- 2w^ sin 20 
\J[a + b — c—d — {a —b—c + d) cos 0]^ -I- 4ic2 sin^ 0 


(42) 


All three possible variants for the quantum discord (Qg, Qii-/ 2 j and Qe) 
can really exist in physical systems. In the case when a = b and b = c (or 
Si = S 2 = 0) the conditional entropy minimum is always achieved at one of 
two bound points [S]. However, this is wrong for the more general X states; 
global minimum can take place at inner points of the interval (0,7r/2). 

Indeed, following the authors eg, let us consider the state 

/ 0.0783 0 0 0 \ 

0 0.125 0.100 0 

0 0.100 0.125 0 

\ 0 0 0 0.6717/ 


(43) 





















On the quantum discord of general X states 


9 



Fig. 2 Quantum conditional entropy Scond a function of measured angle 0 for the state 
I® 


Using Eqs. (I^ - (I5^ we have computed the function Scond{d) for this state. Its 
behavior is shown in Fig. [21 From the figure we conclude that the conditional 
entropy minimum is situated in the intermediate region, namely, at the angle 
9 = 0.4883 ~ 28°. Two other similar numerical examples of quantum states 
are given in Ref. [TH] • 

These examples clearly show that the optimal measurement angles can 
really be in the intermediate region (0,7r/2), i.e., the optimal observables for 
quantum discord can be not only the ax or az, but also their superposition. 

For the real X state with constraint |u + u| > |u — u| (i.e., uv > 0 or 
signu = signu), the authors [17) have proved a theorem which guarantees that 
the optimal observable is az if 

i\u\ + \v\r<{a-b){d-c) (44) 

and ax if 

|u| + kl > iVad— Vbc\. (45) 

The theorem states nothing for the region lying between these bounds. But in 
the case 

ac = bd (46) 

the inequalities (HU) and (HSl) lead to absence of the intermediate region [2S]. 


4 Equations for the boundaries 

Let us start with a heuristic example. Consider a two-parameter family of X 
states [mus] 

/e/2 0 0 e/2\ 

_ 0 (1 — e)m 0 0 

^“0 0 (l-e)(l-TO) 0 

Ve/2 0 0 e/2/ 


(47) 
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P = ^[1 + (1 - e)(2TO - l)(cr? - cri) + e(o-iCr^ - crfcrf) + (2e - l)criCr^] 

= e|<P+)(<l>+| + (1 - e)m|01)(01| + (1 - e)(l - m)|10)(10|, (48) 


where = (|00) + |ll))/-\/2- The given density matrix p represents the 
generalized Horodecki states 

Simple calculation yields Qo = £ (in bits). Sufficient conditions and 
(HH|) for the Qo and < 3^/2 subdomains give [T7] 


^ 2m(l — m) 

^ “ 1 + 2m(l — m) 

and 

^ ■\/to(1 — m) 

1 + Y^m(l — m) ’ 

respectively. But in the region 


(49) 


(50) 


2m(l — m) to(1 — to) 

1 + 2 to (1 - to ) ^ 1 + .^ to (1 - to ) 


(51) 


the above theorem does not say anything. 

Let us now find the lines on the plane (to, e) which are defined by the 
condition 

Qo{m,e) = Q^{m,e). (52) 

Then we will study the changes of curves Scondifi) in the neighborhood to 
those lines. 

Using Eqs. (EZl), (IMl), and (1571) we have numerically solved the transcen¬ 
dental equation (IH51) . The solution is only one. The results are plotted in Fig. [3] 
by dotted line. 

Consider in detail a particular case. Let the e is held fixed and equal, for 
example, to e = 0.228 (see Fig. IS. Then the equality Qo = Q'n /2 is satisfied 
at the crossing point TOx = 0.101234. Study now the behavior of Scond{0) 
when the parameter to varies. Inequalities (1491) and (1501) guarantee that when 
TO < 0.096 545 the discord equals Q = Q ^/2 and Q = Qo when to > 0.180 107. 
If TO = 0.1015, the minimum of Scond{0) is at 0 = 0 [see Fig.|4](a)]. Moreover, 
the angle P = 0 is optimal for all larger values of to. When the to decreases 
the minimum on the curve Scond{d) inside the interval between 0 and 7r/2 
appears. The minimum is clearly seen when to = 0.1014 [Fig. |3](6)]. Near the 
point TO = 0.101 234 the minimum achieves large depth. By further decreasing 
TO, the minimum moves to the bound P = 7r/2 and then it disappears at all. 
Optimal measurements undergo to the angle 9 — -k 12. 

We argue now that both lower and upper boundaries of the interval where 
the optimal angles lie between 0 and 7r/2 are exact, i.e., the intermediate 
minimum of Scond{9) suddenly appears and suddenly disappears. Above all, 
we note that the first derivative of function Scond{9) at 0 = 0 and 7r/2 equals 
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Fig. 3 Subdomains Qr:j 2 -> Qo? (between them) Qq for the state (147D . Dotted line 
corresponds to the condition (5,^/2 = Qo- Solid lines 1 and 2 are the 7r/2- and 0-boundaries, 
respectively 



Fig. 4 Appearance and disappearance of an intermediate minimum on the conditional 
entropy curve by transition from Qq to subdomain. Here, Scondi.^') corresponds to the 

state 1I47I I at the fixed value of e = 0.228 and m = 0.1015 (a), 0.1014 (&), 0.101234 (c), 
0.1011 (d), and 0.1008 (e) 


zero in general case: 5 ''q„^(0) = = 0. This is easy to check by direct 

calculations using Eqs. Let us turn now again to the Fig. 01 By 

fixed value of parameter e and for each value of m one can say at any moment 
the inside minimum exists or it is absent. For instance, when m = 0.1015 
(e = 0.228) the function Scond{0) is concave at the point 0 = 0 and therefore 
its second derivative S”^^^{Tr/2) < 0. But when m = 0.1014 the conditional 
entropy has a local maximum at the same bound point 0 = 0 and therefore 
S"ondi'^/‘^) > 0- Hence, the bifurcation point (in the sense that two extrema 
arise from one ) m is determined by the condition 

ScondiO) = 0. (53) 


Similarly we have for the other bound point 0 = 7r/2, 

Scoudi^/2) = 0. 

Using Eqs. (I^ - (I5T]) we get the second derivatives at limiting points: 


^cond 


(0) = - (a — 6 + c — d) f 2 In —— 

1 / . .M ad 1 2 / 1 

+ -{a -b - c + d)\n- - -w^ - 

4 be 2 \ a — c 


In 


In- 


bd 


1 1 ^ 
737 


(54) 


(55) 
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and 


S’LnA'^m = —[(«- C){b -d) + w^]\n^^+{a-b + c-df 

11 2 

— —T--r \a — h-\-c — d-\—(a -\- h — c — d){a — h — c -\- d)\ 

2(l + r)^ 

11 2 

-7- -\a — h-\-c — d - (ah — c — d){a — h — c-\-d)\ , 

2(1 — r) r ■' 


(56) 


where 

r =[{a + b-c-df (57) 

and w is given by Eq. (15^ . The relations (l55]) - (l571) are the boundary equa¬ 
tions for the crossover subdomain Qg Thus, the boundaries consist of 

bifurcation points. 

If the solutions of Eqs. (15^ and (IMl) are the same, the intermediate subdo¬ 
main Qg is absent and the quantum discord is given by analytical expressions. 
On the other hand, instead of rough conditions (Hill and (HH1) . the inequalities 
S"ondiO) ^ 0 5 '"Qjj^( 7 r/ 2 ) < 0 define now the complete subdomains Qo and 

Q^/ 2 , respectively. 

Numerical solution of Eqs. Q-dnil) for the state dUl) shows that the 
boundaries are the lines going approximately parallel to the dotted lines (see 
the lines 1 and 2 in Fig. [SD- As a result, the subdomain appears within which 
the optimal angles should be found numerically. Out of this subdomain we 
have analytical expressions for the quantum discord. By e = 0.228, the value 
for m of 7 r/ 2 -boundary equals to ^/2 = 0.100 997 and for the 0-boundary it is 
Too = 0.101474. The middle of this interval equals 0.101 236 which is near the 
point mx = 0.101234. 

Consider the discord behavior by a transition from the subdomain Q ^/2 to 
Qo one (Fig. O. One can see that down to crossing point rrix = 0.101 234 the 
discord Q, according to Refs. [aiiniin], equals Q ^/2 and above the point TOx it 
equals Qo (see Fig.[S]). If this were valid, the discord Q = min{Q ^/27 Qo} would 
not be differentiable at the intersection point mx. However, in fact, the true 
discord Q = min{Q.n./ 2 ,Qe,Qo} is smooth. This follows from the numerical 
solution of the task in the intermediate domain. The results are shown again 
in Fig. EJby solid line. It is clearly seen that smoothness occurs. We may say 
that, instead a fracture at mx, two hidden transitions occur at the 7r/2- and 
0 -boundaries. 

Notice that the conditions pg]) and (ISOj) are rough too and lead to the 
bounds which lie far beyond the region of Fig. [31 


5 Bell-diagonal states 

The case a = d and 5 = c or Si = S 2 = 0 (i.e., when both local Bloch vectors 
are zero) corresponds to the Bell-diagonal states. Domain of definition for the 
physical states, 27, lies now in the three-dimensional cube defined by ci, C 2 , C 3 G 
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m 


Fig. 5 Dependencies of the false discord Q = min{( 57 r/ 2 » Qo} (dotted line) and the corrected 
quantum discord Q = min{Q 7 r/ 2 i Q 05 Qo} (solid line), for the state (I47lt with parameter 
e = 0.228 Longer bars mark the exact boundaries m^j 2 = 0.100997 and ttiq = 0.101474. 
Subdomains m < m^j 2 , < mo, and m > mo correspond to the discord branches 

Q-Kj 2 ‘> Qqi Qo respectively 


[— 1 , 1 ]. Two second-order hypersurfaces (EOl) and (EH are transformed to the 
two first-order surfaces 


± Cl -1- C2 -I- C3 - 1 = 0 

(58) 

± Cl - C2 -b C3 -I- 1 = 0. 

(59) 


The former consists of two semi-planes with a A-shaped cross section and the 
latter is similar to it but has a V-shaped cross section. The angle between 
semi-planes equals arccos(l/3) ~ 78°. These semi-plane surfaces put bounds 
to the domain V that is reduced, as shown in Fig. |6l to a tetrahedron with 
vertices [5B] 

Ui = (-1,1,1), U2 = (l,-l,l), U3 = (l,l,-1), U4 = (-l,-l,-l); 

(60) 

these vertices lie in octants II (—,-|-,-|-), IV (-1-,—,-|-), V (-|-,-|-,—), and VII 
(—,—,—), respectively. The centers of tetrahedron facets are 

01 = (1/3,1/3,1/3), 02 = (-1/3, -1/3,1/3), 

03 = (-1/3,1/3,-1/3), 04 = (1/3,-1/3,-1/3). (61) 

Tetrahedron volume equals a third (i.e., about 33.3%) of the cube one. Notice 
that the tetrahedron vertices are the states with maximal value of discord 
(which equals one in bit units). 

It is known [55] that the states with zero discord are negligible in the 
whole Hilbert space. In particular, it has been proved mm that, when si = 
S 2 = 0, the zero-discord states have at most one nonzero component of vector 
(ci,C 2 ,C 3 ), i.e., all classical-only correlated states lie on the Cartesian axes 
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Fig. 6 Tetrahedron with vertices and 'U 4 is the domain for the total Bell-diagonal 

states. Two regions (O, , 1 ^ 2,01 , 02 ) and (O, 'U 3 , D 4 , 03 , 04 ) correspond to the physical states 

with Qo discord 


Oci, Oc 2 or OC 3 . (This corresponds to the so-called “Ising spins” introduced 
as a matter of fact by his adviser W. Lenz in 1920 [5^133) . 1 

In the case of Bell-diagonal states, both boundary equations (ESD-diZl) are 
reduced to a relation 

{a-br = {\u\ + \v\r, (62) 

so that 

2 |c 3 | = |ci-1-C 2 I-I-|ci - C 2 I. (63) 

Thus, the 7 r/ 2 - and 0-boundaries are coincident, the Qg subdomain is absent 
here, and the quantum discord is given by the explicit analytical formula 
Q = minIQoi 0 ^/ 2 } which is in full agreement with Luo’s results [5]. 

From Eq. (1551) . four equations follow 


C 3 = ±ci C 3 = ±C 2 . (64) 

These planes divide the tetrahedron into subdomains Qo and ( 5 ,r/ 2 ! where the 
quantum discord takes the values Qo or QTr/ 2 - Qo subdomain consists of two 
hexahedrons {0,vi,V2,01,02) and {0,V3,V4,03,04); they are shown in Fig. |6l 
The remaining volume of a tetrahedron belongs to the QTr /2 states. It is in two 
times larger than the volume of Qo states. 

The behavior of quantum discord for the Bell-diagonal states along different 
trajectories is illustrated in Fig.[7]by solid lines. Figure[7]a shows the discord as 
a function of C 3 G [—0.95, 0.45] by fixed values of ci = 0.3 and C 2 = 0.25. The 
curve is continuous but has the fractures at C 3 = ±0.3. They happen when the 
trajectory crosses the planes dividing the Qo and Q,r /2 subdomains (see Fig.jS]). 
In this case, the optimal measurement angle 6 varies discontinuously, namely, 
it jumps from 0 = Oto0 = 7 r/ 2 or inversely. In the vicinity of cross points, 
the conditional entropy Scond{d) changes its form going through a straight line 
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Fig. 7 Quantum discord for the Bell-diagonal states: (a), Q = minjQo^ Q7r/2} C 3 by 
Cl = 0.3 and C 2 = 0.25, longer bars mark the positions of fracture points at C 3 = ±0.3; ( 6 ), 
Q = QT^j2 (solid line) and Qq (dotted line) vs C2 when ci = 0.25 and C 3 = 0 


(where any angle 6 G [0, 7 r/ 2 ] is optimal). Such a regime of conditional entropy 
behavior is shown in Fig. [H 

Figure [7] b shows the behavior of branches Qo and Q.^/2 as functions of 
C 2 by fixed values of other two parameters, ci = 0.25 and C 3 = 0. Since here 
Q- 7 r /2 < Co) the quantum discord Q equals Q,r/ 2 - The curve QTr /2 has two 
fractures. This means that the branch ( 57 r /2 is a piecewise-analytic function. 
In this case, however, the optimal measurement angle does not change its value 
0 = 7r/2 and therefore the position of conditional entropy minimum remains 
immutable. 



Fig. 8 Transition between j 2 ^nd Qo subdomains via a strait line for the conditional 
entropy. Here, 5'cond(^) is at and near the fracture point C 3 = 0.3 on the quantum discord 
curve in Fig.The curves 1 , 2 , and 3 correspond to C 3 = 0.29, 0.3, and 0.31, respectively 
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6 Physical systems with the Qg subdomains 

We are interested now in the systems with Qg subdomains. As it was seen 
from the previous section, such regions do not exist in the Bell-diagonal states. 
Therefore in this section we will consider systems with nonzero Bloch vectors. 


6.1 Phase flip channels 

Let us consider the dynamics of quantum discord under decoherence (for a 
recent review, see, e.g, [34] and references therein). The authors [11] have con¬ 
sidered such a dynamics in the phase flip channel. The problem is to calculate 
the quantum discord for the X matrix 

£ = i[l -I- Si(Tz (g) 1 -I- S2I ® CTz + (1 - p)^Ci(Ta, (g) 

+ {1 - p)‘^C2ay (Jy + Csaz ® (Tz], (65) 


Here, the parametrized time p = 1 — exp(—qf), where t is the time and 7 is the 
phase damping rate. The authors m restricted themselves to the case where 

C2 = -C 3 C 1 , S2 = C 3 S 1 , -1 < C 3 < 1, -1 < Si < 1. ( 66 ) 

Expansion coefficients in Eq. (I65|) are related to the corresponding X matrix 
elements as 


a = (1 -I- Si -I- S 2 -I- C 3 )/ 4 , 6 = (1 -p Si - S 2 - C 3 )/ 4 , 

c = (1 - Si -I- S 2 - C 3 )/ 4 , d = (1 - Si - S 2 -b C 3 )/ 4 , (67) 

M = (1 - p)^(ci - C 2 )/ 4 , -y = (1 - p)^(ci -b C 2 )/ 4 . 

Owing to the relation S2 = C3S1, the matrix elements a, b, c, and d satisfy the 
condition (1461) and hence the Qg domain is absent here; conditional entropy 
behaves similar to that as shown in Fig. |5J Thus, nonzero values of si and S 2 
are the necessary but not sufficient condition for existence of Qg phase. 

Consider a different initial state. For example, let us take si = S 2 = 0.65, 
Cl = C 2 = 0.249, and C 3 = 0.5. As seen from Fig. [HI the curves Qoip) and 
Q-kI 2 [v) have a crossing point at px — 3158. An additional study shows that 
the transition Qt^i 2 Qo goes through the appearance of single minimum 
on the Scond{d) curves inside the interval between 0 and 7 r /2 (similarly to the 
curves on Fig. a. 

Solution of equations for the boundaries, Eqs. (ESD-dST]), shows that the 
7 r/ 2 - and 0-boundaries do not coincide now and therefore the Qg region exists 
here (see Fig. uni). 
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p 


Fig. 9 Qo (solid line) and Q-j^ (dotted line) in bits vs p for the phase flip channel with 
parameters si = S 2 = 0.65, ci = C 2 = 0.249, and C 3 = 0.5. Crossing point of the lines is at 
Px = 0.315 789... 



0.3145 0.315 0.3155 0.316 0.3165 0.317 

P 


Fig. 10 Dependencies of the false discord Q = rmn{Q^ 121Q0} (dotted line) and the cor¬ 
rected quantum discord Q = rmn{Q^i2->Q9->Qo} (solid line) vs p for the phase flip channel 
with parameters si = S 2 = 0.65, ci = C 2 = 0.249, and C 3 = 0.5. Longer solid bars mark the 
boundaries Pt ^/2 = 0.314 949 and pQ = 0.316 637. Longer dotted bar marks the position of a 
fracture, px = 0.315 789, on the curve Q{p) 


6.2 Thermal discord 

We now discuss systems at thermal equilibrium. Let us consider the XYZ spin 
Hamiltonian 

H = + Jycrlal + JzCrj'crl + Bial + i?2(T|). 


( 68 ) 
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This Hamiltonian contains five independent parameters J^, Jy, Jz, Bi, B2 G 
(— 00 , 00 ) (i.e., in TZ^) and is the most general real symmetric traceless X 
matrix. The corresponding Gibbs density matrix is given as 


PAB = 


(69) 


(here /3 = 1/T, T is the temperature in energy units, Z is the partition 
function) and has also the five-parameter real X structure. Thus, the map 
{Bi/T,B 2 /T,Jx/T,Jy/T,Jz:/T) o (si, S2, ci, C2, C3) (that is TZ^ O V) allows 
in general to change the density-matrix language on a picture of interactions 
in the XYZ dimer in inhomogeneous fields Bi and B 2 . 

Having solved eigenproblem for the Hamiltonian (1681) we then find expres¬ 
sions for the thermal density matrix elements 


_ 1 cosh(/?i?i/2) -H [(Hi -H B 2 )/Ri] sinh(/3j?i/2) 

2 cosh(/3i?i/2)-|-exp(—/SJj) cosh(/3i?2/2) ’ 

1 cosh(/3i?2/2) -H [(Hi - B 2 )/R 2 ] smh(^j?2/2) 

2 exp(/?Jj;) cosh(/?i?i/2)-|-cosh(/3i?2/2) 

_ 1 cosh(/3i?2/2) - [{Bi - H2)/i?2]smh(^j?2/2) 

2 exp(,dJ2) cosh(/3i?i/2)-I-cosh(/3i?2/2) 

1 cosh(/3j^i/2) - [(Hi -H B 2 )/Ri] sinh(/3i?i/2) 

2 cosh(/?i?i/2)-I-exp(—/3J2) cosh(/3i?2/2) ’ 
_ 1 [{Jx - Jy)/Ri] sinh(/3i?i/2) 

^ 2 cosh(/3i?i/2) -I- exp(—/3J2) cosh(/3i?2/2) ’ 

^ 1 [{Jz, + Jy)/R2]smm3R2/2) 

2 exp(,5Jz) cosh(/3i?i/2) -I- cosh(/3i?2/2) ’ 


where 


Ri = [iB,+B2f + iJ.-Jyf]^/^ R2 = [iBi-B2f + {J. + Jyf]^^^. (71) 
For the correlations functions (USD, we have respectively 

ci^2 = |[±[(J. - Jy)/Ri]e^^^/^smh{l3Ri/2) 

+ [{Jx + Jy)/R 2 ]e-^^‘/^ sinh(/ 3 i? 2 / 2 )], 

C 3 = cosh(/3i?i/2) — cosh(/ 3 i? 2 / 2 )], (72) 

Zj 

si,2 = |[[(Si + 

±[(Hi - H2)/i?2]e-''‘''/2sinh(^i72/2)l, 
where the partition function equals 

Z = cosh(/3i?i/2) -h cosh(/ 3 i? 2 / 2 )] (73) 

and i?i and R 2 are given again by Eq. (ED. 
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Fig. 11 Dependencies of the false discord Q = min{Q 7 r /25 Qo} (dotted line) and the correct 
quantum discord Q = i 2 iQ 9 'iQo} (solid line) for the XXZ dimer with parameters 

J = J 2 = 1.02 and B = 1. Longer bars mark the temperatures T ^/2 = 0.76106 and 
To = 0.85361. Domains T < 71^/2 < < ^0 ) 7" > Tq correspond to the discord 

branches Qt^i 2 ^ Qq, s,nd Qq, respectively 


For each choice of interaction constants Jx^ Jy^ Jz cind external fields Bi 
and B 2 we will find the points where the condition Qq = Qtt /2 is satisfied. After 
this we will again study the changes of curves Scond{0) in the neighborhood of 
points found. 

Taking, for example, a dimer with parameters Jx = Jy = J = Jz = 1 - 02 , 
and Bi = B 2 = B = 1 (that is the XXZ dimer in an uniform field) we consider 
the thermal discord behavior by a transition from the subdomain QTr /2 to Qo 
one ('Fig. fTTl) . From the figure one can see that down to the crossing point Tx = 
0.81296 the discord Q, according to Refs. [aiiniin], equals QTr /2 and above 
the point Tx it equals Qo- If this were valid, the discord Q = miTi{Qo,Q^/ 2 } 
would have a fracture at the intersection point Tx. However, in fact, the true 
discord Q is a smooth function (at least, it is a function of differentiability class 
C^). This follows from the numerical solution of the task in the intermediate 
domain, where the Scond(Q) curves change similar as in Fig. |4l Results for 
the quantum discord are shown again in Fig. 1111 At the bifurcations points 
Tt ^/2 = 0.76106 and Tq = 0.85361, the higher derivatives of quantum discord 
Q — min{Q,n./ 2 , Qe, Qo} exhibit a discontinuous behavior. 

6.3 Heteronuclear systems with dipolar coupling 

Let us consider the system (1^ with parameters Jx = Jy = —D and Jz = 2D. 
Such a model corresponds to a dipolar coupled dimer which is stretched along 
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BijD 


Fig. 12 Subdomains Qt^/ 2 -) Qo-, and (between the lines 1,2 and 1^2') Qq for the spin 
dimer 111 at the normalized temperature T/D = 1. Dotted lines 3 and 4 correspond to 
B 2 = 4Bi and B 2 = BijA.^ respectively 


the z axis [3S] 


n = + a\al - 

Here the dipolar coupling constant (in frequency units) equals 

„ ^ /ip 7172 

47r 2rQ 


(74) 


(75) 


where is the magnetic permeability of free space, 71 and 72 are the gyro- 
magnetic ratios of particles in the dimer, and tq is the distance between those 
particles. Normalized fields Bi and B 2 in Eq. (ESI) are 

Bi = 7 iBo, B2 = 72.B0, ( 76 ) 

where Bq is the external magnetic field induction. 

We have performed necessary calculations (according to our approach de¬ 
veloped in the previous sections) and found the subdomains of quantum dis¬ 
cord in the plane {Bi/D,B 2 /D). The results are shown at the normalized 
temperature T/D = 1 in Fig. [121 From this figure one can see that such a sys¬ 
tem has the Qe regions (between the 1,2 and 1', 2' lines). These regions can be 
reached by varying the external magnetic field Bq. Two possible trajectories 
are shown in Fig. HU by dotted lines, B 2 = 4i?i and B 2 = i?i/4. (The value 
72/71 = 4 approximately corresponds to the quotient of gyromagnetic ratios 
for the nucleus of and ^^C.) 
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We found also that in the Qg subdomain the conditional entropy Scond{0) 
has only one minimum that is located in the interval (0,7r/2). The picture is 
qualitatively similar to that is shown in Fig. |4l 

So, the Qg region and corresponding sudden changes of quantum corre¬ 
lation behavior at their boundaries can be observed in solid materials with 
nuclear dimers. 


7 Results and perspectives 

In light of the above, the calculation of quantum discord of any X states can 
be achieved by following steps. First, the density matrix o is transformed to 
the real form, i.e., the quantities u and v are calculated using Eqs. © and 
d?]). It is also well to solve the equation Qq — Qt^I 2 and determine possible 
crossing points of branches Qo and Q^/ 2 - Then the equations = 0 and 

*^cond('^/2) = 0 are solved to find the boundaries for the intermediate subdo¬ 
main Qg. After this, one should numerically find the optimal measurement 
angle 9 G (0,7r/2) and compute Qg = Q{9). As a result, the quantum discord 
is given by Q = min{Qo, Qe, Qt,/ 2 }- 

The formula for calculating the quantum discord belongs to a piecewise- 
defined type 

{ F{x,a), X e fia 

F{x,b), xeHb (77) 

minag(a_t,) F{x, a), x G I2c- 

In other words, the domain of definition, 12, of the function f{x) consists of 
subdomains in which the function is given by closed analytical expressions or 
it exists in a numerical form. 

So, the quantum discord of X states is represented analytically if the Qg 
subdomain is absent. Then the quantum discord is given by the closed form 
Q = min{(3,n./2, Qo}- The discord is continuous, but generally speaking it is 
a piecewise smooth function. In particular, this is valid for a spacial class of 
X states, namely, for the Bell-diagonal states. For them, we found the Qo 
and QTr /2 regions in the total domain of their definition (Fig. 15]). It would be 
interesting to find the subdomains Qo: Qe, and QTr /2 in the five-dimensional 
domain 22 making, e.g., an atlas of maps. 

Also, we have shown in this paper that the boundaries for the transition 
subdomain from Qq to Q 7 r /2 or reversely are exactly defined. They consist of 
nonanalyticity points which are bifurcation ones. The corresponding equations 
for these boundaries have been found. The boundaries may coincide and then 
the quantum discord is evaluated analytically in the total domain of dehnition. 
The regions Qg with the optimal intermediate angles 9 G (0,7r/2) have been 
found for a number of physical systems including the phase flip channels, spin 
dimers at the thermal equilibrium, heteronuclear systems with dipolar inter¬ 
action. The transitions Q ^/2 Qg Qo occur continuously and smoothly. 
This is a new type of transitions for the quantum discord. 
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We have found only two regimes for the conditional entropy change by 
above transitions: (i) via the birth of one intermediate minimum (as shown in 
Fig. |4]) and (ii) via the strait line (as shown in Fig. |8]). It is hoped that our 
observations will be rigorously proofed and, maybe, generalized in the future. 

At present the attempts are made to obtain analytical formulas for the 
super quantum discord of X states with nonzero Bloch vectors [551157] . In 
this connection one should note the the authors do not take into account a 
possibility of intermediate optimal angles for the weak measurements which 
are a generalization of projective ones. 
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